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Background: The available data for E2 transition strengths in the region between neutron deficient hafnium 
and platinum isotopes is far from complete. More and precise data are needed to enhance the picture of structure 
evolution in this region and to test state-of-the-art nuclear models. In a simple model, the maximum collectivity 
is expected at the middle of the major shell. However, for actual nuclei particularly in heavy mass regions, which 
should be highly complex, this picture may no longer be the case, and one should use a more realistic nuclear 
structure model. We address this point by studying the spectroscopy of Hf as a representative case. 

Purpose: To re-measure the 2+ half-lives of 1 ' 2 > x ' 4 - 176 Hf, for which there is some disagreement in the literature. 
The main goal is to measure, for the first time, the half-lives of higher lying states of the rotational band. The 
new results are compared to a theoretical calculation for absolute transition strengths. 

Method: The half-lives were measured using 7-7 and conversion electron -7 delayed coincidences with the fast 
timing method. For the determination of half-lives in the picosecond region the generalized centroid difference 
method was applied. For the theoretical calculation of the spectroscopic properties, the interacting boson model 
is employed, whose Hamiltonian is determined based on microscopic energy-density functional calculations. 

Results: The measured 2± half-lives disagree with results from earlier 7-7 fast timing measurements, but are in 
agreement with data from Coulomb excitation experiments and other methods. Half-lives of the 4 4 and 6 + states 
were measured, as well as a lower limit for the 8 ^ states. 

Conclusions: This work shows the importance of a mass dependent effective boson charge in the interacting 
boson model for the description of E2 transition rates in chains of nuclei. It encourages further studies of the 
microscopic origin of this mass dependence. New experimental values on transition rates in nuclei from neighboring 
isotopic chains could support these studies. 

I. INTRODUCTION 

The absolute strength of E2 transitions between low- 
lying states of even-even nuclei is an important observ¬ 
able to test nuclear models that describe collective phe¬ 
nomena. Well deformed even-even nuclei exhibit large 
quadrupole transition strengths B(E2;2+ —» Of) PQ. In 
a simple picture this strength increases smoothly for in¬ 
creasing number of valence nucleons or holes along an iso¬ 
topic or isotonic chain as collectivity increases. Assuming 
symmetry of particles and holes this model yields a max¬ 
imum at mid-shell [5] . Microscopic effects can break the 
particle-hole symmetry and lead to a different picture. 

Recent measurements on tungsten and osmium iso¬ 
topes EH5] showed that for a given isotopic chain the 
maximum of the absolute strength B(E2; 2^ —is 
not found at mid-shell but at lower neutron number. 

Furthermore there seems to be a sudden increase at N 
= 98, which is not expected for a collective observable. 

The overall picture of available data on B(E2)s is rather 
erratic for these and neighboring nuclei, which makes 
a comparison to nuclear structure models difficult. In 
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many cases important data are missing or there exist dis¬ 
agreeing results from experiments using different meth¬ 
ods. For a better understanding of nuclear structure 
evolution in this region it is important to have a more 
complete picture of absolute transition strengths for the 
isotopes with 70 < Z < 78. 

In this paper, new B{ E2) values for hafnium isotopes 
(Z=72) are presented. They were measured using the 
method of delayed coincidences with LaBr 3 (Ce) detectors 
and an Orange conversion electron spectrometer. The 
results improve the data situation and enable us to test 
current model predictions of both absolute and relative 
transition strengths. 

On the theoretical side, a quantitative and detailed 
description of spectroscopic properties for heavy nuclei 
has been provided by fully-microscopic many-body the¬ 
ories, which include the large-scale shell-model [6H8] and 
the self-consistent mean-field mm, or the energy den¬ 
sity functional (EDF), approaches with a suitably chosen 
effective nucleon-nucleon interaction. 

In the case of strongly deformed heavy nuclei, like the 
ones considered here, the dimension of the shell-model 
configuration space becomes exceedingly large requiring 
an appropriate truncation scheme to reduce the compu¬ 
tational cost while keeping the essential features of low- 
lying collective states. Along this line, one could assume 
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that the nuclear many-body problem is approximated by 
a system of interacting bosons: the interacting boson 
model (IBM) [T2I. The essence of the model is the as¬ 
sociation of the collective nucleon pairs relevant for the 
low-energy states, e.g., with spin J = 0 + ( S ), 2 + (D), 
..., to the equivalent bosonic degrees of freedom (s, d, 
...) by means of a well defined mapping procedure. The 
physical states and their decay properties are obtained 
from the calculations in the boson space dam]. 

In recent years, the IBM has been used with consider¬ 
able success in spectroscopic studies [IMS] of a large set 
of nuclei with proton and neutron numbers in the range 
Z = 50 — 82 and N = 82 — 126, respectively. However, 
those studies are rather phenomenological because the 
parameters of the IBM Hamiltonian are determined by a 
fit to known experimental data. 

On the other hand, the mean-field framework has also 
been successfully applied to the study of nuclear struc¬ 
ture phenomena (9]. Both in the non-relativistic [20] [21] 
and relativistic da nn regimes, it is possible to ob¬ 
tain a global and reasonable description of the ground- 
state properties and collective excitations of all nuclei 
across the nuclear landscape with a single parametriza- 
tion of the corresponding EDF. However, to describe 
in detail spectroscopic properties one needs to go be¬ 
yond the mean-field level to restore symmetries broken 
in the mean field approximation as well as to take into 
account fluctuations with collective coordinates. Much 
effort has been devoted to increase the feasibility of such 
calculations, mostly in the framework of the pure gener¬ 
ator coordinate method [9j l23l - l26j or its approximations 
EZH 2 B]- The results confirm the usefulness and reliability 
of these EDF-based approaches for the study of nuclear 
spectroscopy. 

To investigate the spectroscopic properties of the 
neutron-deficient hafnium isotopes, we use the procedure 
of Ref. [29], that determines the Hamiltonian of the IBM 
from EDF-based mean-field results. The idea behind this 
method is to map the deformation energy surface result¬ 
ing from a set of the constrained mean field calculations 
onto the equivalent energy surface for the system of in¬ 
teracting bosons, that is, onto the expectation value of 
the IBM Hamiltonian in the boson condensate state m- 
The parameters of the IBM Hamiltonian determined by 
this procedure do not require any additional adjustment 
to experimental data. The resulting Hamiltonian is used 
to calculate energy spectra and transition rates. So far, 
the predictive power of the method has been verified in 
all the possible regimes of low-energy quadrupole collec¬ 
tive states: spherical vibrational [E], 7 -soft [33133] and 
well-deformed rotational [33] nuclei. 

Given the predictive power of our procedure it is in¬ 
teresting to check whether the IBM Hamiltonian deter¬ 
mined from the microscopic mean-field calculation can 
explain the absolute transition rates for the neutron de¬ 
ficient hafnium isotopes. The microscopic input used 
is the Gogny energy density functional with the recent 
parametrization DIM [33] , which has been shown [33- 


133 to have a similar level of predictive power in the 
description of nuclear structure phenomena as the stan¬ 
dard parametrization D1S |38l . However, and to con¬ 
firm the robustness of our results, we also discuss the 
results obtained with the more traditional Gogny D1S 
parametrization. 

This paper is organized as follows. In Sec. [TT] we de¬ 
scribe the experimental procedure and the data analy¬ 
sis. We then present the results of the experiment in 
Sec. |III| In Sec. [IVj we compare the present experimen¬ 
tal results with the predictions by the IBM-2 calculation 
combined with the self-consistent mean-field method us¬ 
ing the Gogny energy-density functional. A conclusion is 
given in Sec. [V] Finally, in Appendix [A] the theoretical 
procedure is described in detail. 


II. EXPERIMENTAL PROCEDURE AND DATA 
ANALYSIS 

The experiments were performed at the Institut fur 
Kernphysik of the University of Cologne. The Cologne 
FN tandem accelerator delivered an a beam to induce 
the reactions 170 Yb(a, 2n) 172 Hf and 172 Yb(a, 2n) 174 Hf 
at a beam energy of 27 MeV, and 174 Yb(a, 2n) 176 Hf at 
a beam energy of 26 MeV. The target thickness was 0.40 
mg/cnr 2 , 0.40 mg/cm 2 , and 0.42 nrg/cm 2 , respectively. 
Thin targets were chosen in order to minimize energy 
straggling and absorption of emitted internal conversion 
electrons (ce). These were measured with the Cologne 
Orange spectrometer, a toroidal magnetic spectrometer. 
The electrons are measured with a fast plastic scintillator 
detector at the exit slit of the spectrometer. See refer¬ 
ence [39] for more details on this instrument. A scan with 
different magnetic fields, i.e. different electric currents, 
yields an electron momentum spectrum with the different 
conversion lines corresponding to each nuclear transition. 
Figure [T] shows ce spectra measured for the three inves¬ 
tigated reactions. The exponential background at low 
energies is due to 6 electrons, which are produced by the 
beam ions traversing the target. Note that these elec¬ 
trons are not correlated in time with the decay radiation 
emitted after a nuclear reaction. 

The 7 -rays were measured using a small array of six 
LaBr 3 (Ce) (hereafter called LaBr) scintillator detectors 
and one high purity germanium detector (HPGe) (see 
Fig. § which were mounted perpendicular to the beam 
next to the target position of the Orange spectrometer. 
The LaBr-crystals were cylindrical and 1.5” x 1.5” in 
size. Four of the LaBr detectors were equipped with 
bismuth germanate (BGO) Compton shields and coni¬ 
cal lead collimators to provide active Compton suppres¬ 
sion and passive shielding. The reduction of Compton 
events generated in the LaBr crystal and the active sup¬ 
pression of stray 7 rays produced by primary Compton 
scattering in the experimental surrounding is important 
in delayed coincidence timing measurements. Such back¬ 
ground events, including inter-detector Compton scat- 
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FIG. 1. Conversion electron momentum spectra of the three 
investigated nuclei (7 is the current applied to the magnetic 
spectrometer, see text for more details). Different conversion 
lines from the same nuclear transition are marked accordingly. 
The step intervals during the scans were 0.3 A ( 1,J Hf and 
174 Hf) and 0.5 A ( 176 Hf). 


tering (cross-talk events), are time correlated and the 
determination of their timing contribution poses a major 
source of uncertainty of the final result. The two remain¬ 
ing LaBr detectors were unshielded. The HPGe detector 
was installed perpendicular to the beam axis. Its main 
purpose during the experiment was to monitor the reac¬ 
tion, taking advantage of its energy resolution, which is 
superior to that of the LaBr detectors (see Figure |3](a)). 
In the analysis the HPGe data was used to confirm the 
level schemes and to identify contaminations which could 
otherwise possibly be overlooked in the LaBr spectra. 
This was especially important for the 172 Hf data, as 172 Hf 
is /3 unstable, albeit with a long half-life of 1.87 years [40] . 

Lifetimes were deduced with the method of delayed 
coincidences. The time difference between two signals 
was measured using time to amplitude converters (TACs) 
arranged in a fast timing circuit as described in [41] . The 
data were recorded triggerless in a list mode format and 
analyzed offline. This way it was possible to sort the 
data with a double (cey and 77 ) as well as with a triple 
coincidence condition (ceyy). LaBr 77 projection spectra 
for both cases are shown in Fig. [3](b,c). 



FIG. 2. (Color online) Photograph of the 7 detector array in 
maintenance position. The position of the four BGO shielded 
LaBr detectors (red circles) and the two unshielded LaBr de¬ 
tectors (blue crosses) can be seen. Only the BGO shield is 
mounted where the germanium detector was positioned dur¬ 
ing the experiment (green star). The view is from the target 
position, the beam direction is from right to left. 


A. Half-life of 2+ states 

The half-life of the 2f state of 172 Hf and 174 Hf was 
determined using ce -7 coincidences. For this purpose 
the Orange spectrometer was set to the electric current 
corresponding to the L-conversion line of the 2f —» Of 
transition. The L-conversion line was preferred over 
the K-conversion line because the latter lies at 23 keV, 
where it is buried in the (5-electron background, which 
increases exponentially towards lower energies. Further¬ 
more, the K-conversion coefficient «k is smaller than 
a l for the 2f -A Of transition in the investigated nu¬ 
clei. E.g. a = 5.77, ok = 1.18, ql = 3.49 for an 88.35 
keV E2 transition in 176 Hf. A LaBr energy gate was 
then set on the 4+ -A 2f transition to produce the TAC 
spectra shown in Figure H (a,b). Only the Compton sup¬ 
pressed LaBr detectors were used. A fit of the slope yields 
the half life of the first excited 2 + state. Several differ¬ 
ent parameterizations of the random background were 
tried. All variations were consistent within the uncer¬ 
tainty. The results are shown in Table 0-nm 

In the 176 Hf experiment we encountered a problem 
with the TAC that was started by the electron detec¬ 
tor which was not noticed until after the beam time. It 
was not possible to extract reliable ce -7 time spectra and 
the 7-7 approach was therefore used in this case. The 
2+ -A Of transition is highly converted and at very low 
energy and the 7 line in the LaBr spectra is therefore 
weak and on top of a lot of time-correlated background. 
The background subtraction, which yields the spectrum 
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FIG. 3. 7 -ray spectra of 174 Hf from one experimental run 
with different coincidence conditions, (a) A spectrum mea¬ 
sured with the Ge detector in coincidence with electrons corre¬ 
sponding to the L-conversion of the 2+ —► 0+ transition. The 
Ge detector was shielded with lead and copper to suppress X 
rays, (b) LaBr 7-7 projection without electron coincidence 
condition, (c) LaBr 7-7 projection with the same electron 
coincidence condition as in (a). 


shown in Figure[4](c), increases the uncertainty dramati¬ 
cally. The result agrees well with the value from Coulomb 
excitation measurements given in the nuclear data sheets 
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FIG. 4. ce -7 (a,b) and 7-7 (c) delayed coincidence time spec¬ 
tra of the cascade 4+ —> 2+ —>• 0J~. Compton background has 
been subtracted to yield spectrum (c). The half-life was de¬ 
termined by fitting the slope. The best fit of an exponential 
function is shown as a solid line. The random background, 
shown as a dashed line, was determined separately. 


B. Half-life of higher lying yrast states 

For the measurement of the, previously unknown, 
shorter half-lives of the 4 + , 6 + , and 8 + states, TAC spec¬ 
tra of 7-7 coincidences between two LaBr detectors were 
analyzed. LaBr coincidence energy spectra are shown in 
Figure [5] For these fast timing measurements in the ps 
region the generalized centroid difference (GCD) method 
was employed 0D> a refinement of the centroid shift 
method [43]. With the GCD method, the centroids of 
two time distributions are measured for a combination 
of transitions that mark the population (feeding transi¬ 
tion) and the depopulation (decay transition) of a nu¬ 
clear excited state. In the case of the present setup, 
if the half life of the excited state is shorter than lps, 
the decay can be considered prompt. If the life time is 


longer a shift of the centroid is observed. If the decay 
is gated on the stop (stop signal on the TAC), the re¬ 
sult is a delayed time spectrum with a centroid Cd e i a yed 
shifted to the right with respect to the prompt distri¬ 
bution. If the decay is gated on the start (start signal 
on the TAC), the result is an antidelayed time spectrum 
with centroid Cantideiayed- A lifetime will show itself as a 
shift of Cantideiayed to the left with respect to the prompt 
position. The halflife of the excited state can be mea¬ 
sured by measuring the difference between the delayed 
and antidelayed centroid if the centroid difference of two 
prompt signals for the respective energy combination of 
decay and feeder is known. The necessary calibration 
of the prompt response difference (PRD) was obtained 
from measurements with a 152 Eu source, as described 
in detail in [5T], using the calibration function given in 
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FIG. 5. Examples of gated LaBr 77 coincidence spectra for 
174 Hf. The transition which was gated on is indicated in the 
top-right corner of each panel. No coincident electron was 
demanded. 


the same reference. Figure [ 6 ] shows the PRD calibration 
points and the fitted function for the reference energy 
344 keV for the two parts of the experiment. An ac¬ 
curacy of ± 10 ps was adopted. Given good statistics 
and peak-to-background ratio, this is the main contri¬ 
bution to the experimental uncertainty. The experimen¬ 
tal observable to be measured is the centroid difference 
AC = Cdeiayed-Gantideiayed which is determined from the 
delayed and antidelayed time spectra like the ones shown 
in Figure [7] The difference between AC and the mean 
PRD, shifted to a given reference energy, corresponds to 
twice the lifetime r = jAS-. 

2 r = |AC-PRD(AC 7 )| (1) 

PRD(A£’ 7 ) = PRD (Cf ee der — -Edecay) (2) 

Analyzed time spectra which are measured using 7 de¬ 
tectors, like the ones shown in Figure [7] necessarily con¬ 
tain background from Compton scattering and crosstalk 
events. These are correlated in time with the transitions 
gated on, but not in the same way as the full energy 
peaks. This timing background is delayed with respect to 
the prompt response for full energy events which leads to 


a systematic error in AC. Therefore a time correction re¬ 
lated to the Compton background contribution was per¬ 
formed in all cases, according to the procedure outlined in 
reference [23]. The centroid of the background distribu¬ 
tion at the gate position was determined by interpolation. 
By assuming that the measured distribution is the sum of 
the background distribution and the desired full energy 
response, the centroid of the latter can be determined 
with the additional knowledge of the peak-to-background 
ratio (ptb). The correction of course contributes to the 
measurement uncertainty or 1/2 . That is 


°r 1/2 



<t(AC- ACq ,) ] 2 

ptb 


ln( 2 ) 

2 


(3) 


with a ac from the measurement of the centroid differ¬ 
ence, and crpRD = lOps from the PRD calibration, and 
the uncertainty from the background correction, which 
includes the interpolated centroid difference of the Comp¬ 
ton background ACc.- The uncertainties are those for 
the actually measured quantity 2r. The result is then 
scaled to the half-life value via the rightmost factor in 
Eq. [3] Dependent on the ptb and the time shift between 
the measured CD and the interpolated time response of 
the Compton background ACc., the correction can be 
important. For an estimation of the magnitude of the 
correction and its effect on the uncertainty of the final 
half life we give an example in numbers: In the case of 
the 311-401 coincidence in 1,4 Hf ptb — ratio ~ 5. The 
absolute value of the correction of AC is 2 ps, i.e. rather 
small. The uncertainty of T 1 / 2 is <tt 1/2 = 3.9ps without 
background correction and <tt 1/2 bg = 4.8ps with back¬ 
ground correction. In the case of the 311-417 keV coinci¬ 
dence in 174 Hf, with ptb — ratio ss 2 , the absolute value 
of the correction of AC is about 16 ps. The uncertain¬ 
ties are a t 1/2 = 4.Ops and ot 1/2 bg = 5.9ps. From this 
example it can be seen, that the background correction 
is very important, especially in the second case. The ef¬ 
fect on the uncertainty of the final half life, however, is 
small compared to the contribution of the overall PRD 
uncertainty of 10 ps. 

It was possible to set a coincidence condition on the de¬ 
tection of a 2 + —► 0 ^ electron for the 7-7 timing analysis 
to clean the LaBr spectra. This improved the peak-to- 
background ratio in the LaBr energy gates used for the 
GCD analysis by a factor of up to 1.4 in the energy region 
below 300 keV. An example is shown in Figure [3] The 
statistics are of course drastically reduced. In the case 
of the strong 4 + —► 2 + -transition, however, the system¬ 
atic uncertainty due to the background is the most severe 
among all transitions. Therefore this approach was only 
used for the determination of the 4+ half-life. Through¬ 
out the analysis only time spectra from coincidences be¬ 
tween shielded LaBr detectors were used, except in the 
determination of the lower half-life limit of the 8 ^~ states. 
Here the transition energies are greater than 400 keV, 
i.e. above the region where crosstalk events play a sig¬ 
nificant role [431 . The measurement of the centroid of a 
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FIG. 6. PRD calibration points measured with a 152 Eu source 
for m Hf, 174 Hf (a), and 176 Hf (b). The reference energy is 
E 1 = 344 keV. The fitted calibration function is shown as a 
solid line. The uncertainty of ±10 ps is indicated as shifted 
dashed lines. For better comparison the data point at E 7 = 40 
keV is shown as an inlay. See text for more details. 


TABLE I. Measured half-lives in 172 Hf. If there is more than 


one intermediate level, and where the final result is an upper 


limit, the effective half-life is given in the column Th/2- The 
adopted half-life values are set in bold face. The values de¬ 
termined with the GCD method are corrected for background 
contributions. For more details see text. 


state 

Edecay 

(keV) 

Efeeder 

(keV) 

T’l/a 

(ps) 

T 1/2 

(ps) 

2 7 

95 

214 


1250(40) 

4 7 

214 

319 


66(5) 


214 

409 

78(6) 

63(12) 


214 

484 

88(13) 

70(17) 


214 

543 

95(18) 

73(20) 

67 

319 

409 


15(8) 


319 

484 

18(9) 

15(12) 

87 

409 

484 

3(7) 

< 10 



TABLE II. Same as in 

Table 0 but for 

■ 174 Hf. 


Edecay 

Efeeder 

TT/a 

Tl/2 

state 

(keV) 

(keV) 

(ps) 

(ps) 

2 1 

91 



1280(40) 

4 1 

207 

311 


77(5) 


207 

401 

101 (6) 

85(9) 


207 

476 

99(15) 

84(15) 


207 

941 

65(47) 

49(52) 


207 

1252 


72(47) 

67 

311 

401 


16(5) 


311 

476 

22 (6) 

18(8) 

87 

401 

476 

5(5) 

< 10 


single distribution imposed an uncertainty between ± 0.5 
ps and ± 4 ps, depending on available statistics. Time 
spectra of direct feeder-decay cascades for the 4+, 6^", 
and 8+ states are shown in Figure [ 7 ] The spectra were 
re-binned to 8 ps per channel for the purpose of display. 
The original resolution was 2 ps per channel. All mea¬ 
sured half-lives are summarized in the Tables m-un 

One merit of the GCD, which it inherits from the cen¬ 
troid shift method [431, 46], is, that lifetimes ti,T 2 ,... 
of intermediate levels connected by a cascade of several 
transitions lead to a centroid shift about an effective life¬ 
time t' which corresponds to the sum of the single life¬ 
times 05] 


t' = t 1 +t 2 -\ -. 


This way it is possible to check the measured half-lives, 
and not least the PRD calibration, for consistency by an¬ 
alyzing time spectra with several coincidence conditions, 
i.e decay-feeder combinations. Where possible this ap¬ 
proach was also used (see Tab. a-®- 


III. EXPERIMENTAL RESULTS 


All measured half-lives are shown in Tables |IV| along 
with B(E2) values. Conversion coefficients for the deter¬ 
mination of B(E2) values were calculated using BrlccFO 
07], see Tab. [v] 

Figure [8] shows a comparison of the results from 


TABLE III. Same as in Tabic 0 but for 176 Hf. 


state 

Edecay 

(keV) 

Efeeder 

(keV) 

T’l/a 

(ps) 

Tl/2 

(ps) 

27 

88 

202 


1470(60) 

47 

202 

307 


90(6) 


202 

401 

110 (6) 

93(9) 


202 

483 

114(6) 

92(12) 


202 

736 

109(8) 

93(11) 


202 

1043 


97(13) 

67 

307 

401 


17(6) 


307 

483 

23(7) 

18(11) 

87 

401 

483 

6(9) 

< 15 
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FIG. 7. (Color online) 7-7 time difference spectra for direct decay-feeder cascades in the nuclei ll2 Hf (left), 1,4 Hf (center), and 
1,6 Hf (right). The delayed spectra, with the decay gated on the stop branch, are shown in gray (color: green). The antidelayed 
spectra, where the decay is gated on the start branch, are shown in black (color: red). The top row shows the time spectra for 
the 6 ) 1 " —> 4+ —) 2+ cascade (with ce gate), the center row those for the 8 ) 1 " —> 6 ^ —> cascade, and the bottom those for the 
lOj 1 " —> 8 ^ —> 6 ^ cascade (both without ce coincidence condition). 




240 

176 H f 


220 


Clx 

200 

. . 

P-ce T Clx 

--T - - cix- ------ -y- -R atocn - - 

180 

1 

I 1 

160 



140 

120 




this Fos Ron Ham Han Tan 

work 1963 1977 1973 1961 1984 


FIG. 8 . (Color online) Comparison of measured B(E2,2~t —> 0+) values of 1 T 2 Hf, 174 Hf, and 176 Hf. The respective adopted 
values found in the nuclear data sheets [401 1421 [481 are shown as a solid line with dashed lines indicating the uncertainty. Other 
data are from Hanl961 [49], Bjel963 [50], Fos 1963 [51], Abol967 [52], Chal971 [S3, Ejil971 [54], Ronl977 [55], Haml973 
[55] . Tanl984 [57] . and Reg2011 03- The experimental method is indicated by ce-ce, 7 - 7 , (3-ce (delayed ce-ce, 7-7 and /?-ce 
coincidence fast timing), Clx (Coulomb excitation), and /r-atom (muonic atom spectroscopy). 

























































TABLE IV. Adopted values of the measured half-lives (see Tables and the corresponding B(E2) value for the transition 
to the next lower lying state in the ground state rotational band. 




m Hf 


m H f 


176 Hf 


Tl/2 

B(E2) 

Tl/2 

B(E2) 

Tl/2 

B(E2) 

state 

(ps) 

(W.u.) 

(ps) 

(W.u.) 

(ps) 

(W.u.) 

2 t 

1250(40) 

194(6) 

1280(40) 

199(6) 

1470(60) 

182(7) 


66(5) 

274(18) 

77(5) 

270(20) 

90(6) 

251(18) 

6 1 

15(8) 

177(1S) 

16(5) 

197(1®°) 

17(6) 

195(tg° 6 ) 

8 1 

<10 

>85 

<10 

>91 

<15 

>60 


this work with results from earlier measurements of 
B(E 2,2+ —>• Of) in the three hafnium isotopes which 
were investigated. For 172 Hf there is one other measure¬ 
ment which also applied the 77 fast timing method, by 


Abou-Leila m who gives a value of T 1 / 2 = 1.55(10) ns. 
This results in B(E2 , 2+ —► 0+) = 154(11) W.u., which is 
several crs lower than the value from this work. Another 
measurement, performed using ce-ce fast timing with 
the Cologne Double Orange Spectrometer, resulted in 
B(E2,2+ —$■ Of) = 199(3) W.u. 03], in agreement with 
the value found in this work. The situation is similar for 
1 , 4 Hf, for which the Nuclear Data Base gives an adopted 
value of B(E2,2+ -> Of) = 152(8) W.u. gS] . Two 
sources are cited which also used 7-7 fast timing. From 
those the adopted value is calculated. One measurement, 
however, used Coulomb excitation [S3], yielding a sig¬ 
nificantly larger value of B(E2,2+ —► Of) = 184(13) 
W.u. This is in agreement with the value determined in 
this work. A ce-ce fast timing measurement resulted in 
B(E2,2f —» Of) = 201(3) W.u., also much higher than 
the value adopted by the Nuclear Data Sheets. For 176 Hf 
the value of this work agrees well with the one adopted 
by the NNDC. In this case the sources do not contain any 
7-7 fast timing experiments from before 1972. The mea¬ 
surements which resulted in the low values cited above 
were all performed in the late 1960s, and used Nal de¬ 
tectors. The energy resolution of these detectors is much 
worse than that of LaBr detectors used today. Under 
such circumstances it is difficult to control and estimate 
the background contributions, especially for transitions 
as low in energy as the 2f —> Of transition in deformed 
nuclei, which typically lie around or below 100 keV. For 
this reason we will adopt the B(E 2, 2f — > Of) values 
measured in this work for our further discussion. The 
other half-lives and absolute transition strengths mea¬ 
sured in this work are measured for the first time. 

The new B{E2) values complete the picture of evo¬ 
lution of quadrupole transition strength in the ground 
state band of mid-shell hafnium isotopes. The result 
is a smooth increase from low neutron numbers N to¬ 
wards mid-shell and a subsequent smooth decrease (see 
Fig. (a)). It is well established that the inves¬ 
tigated hafnium isotopes display the characteristics of 
an axially deformed rotor. In such nuclei the ratio 
-B 42 = B(U2,4f —► 2f)/B(E2, 2f —> Of) can be calcu- 


TABLE V. Internal conversion coefficients (ICCs) from 
BrlccFO 07] used for the calculation of B(E2) values. 



internal 

conversion coefficient a 

transition 

172 Hf 

174 Hf 

176 Hf 

21 ^ 0 + 

4.32(6) 

5.12(8) 

5.86(9) 

4+^2+ 

0.230(4) 

0.256(4) 

0.278(4) 

6+^4+ 

0.0660(10) 

0.0711(10) 

0.0739(11) 

8 +^ 6 + 

0.0327(5) 

0.0345(5) 

0.0345(5) 


lated as the ratio of Clebscli-Gordan coefficients (Alaga 
rules), which yields R 42 = 1.43. As can be seen in Fig- 
ure 0 ](b), the newly measured transition strengths fit well 
into the systematics for rigid rotors, as expected. 


IV. CALCULATION AND DISCUSSION 

A. Description of the model and discussion of 
mean-field results 

To help understand the data, we have performed the¬ 
oretical calculations of spectroscopic properties by em¬ 
ploying the recently proposed methodology of Refs. {2]J, 
[33] . The essential idea of the method is to determine 
the Hamiltonian of an appropriate version of the IBM 
by computing the bosonic deformation energy surface so 
that it reproduces, in a way described below, the basic 
topology of the deformation energy surface of the many- 
fermion system computed microscopically with the HFB 
method. The resultant Hamiltonian is used to calculate 
energy levels and wave functions of excited states. 

Our starting point is the microscopic calculation of 
the deformation energy surface within the self-consistent 
mean-field model. We have performed, for each individ¬ 
ual Hf nucleus, a set of the constrained Hartree-Fock- 
Bogoliubov (HFB) calculations, and obtained the defor¬ 
mation energy surfaces in terms of the quadrupole col¬ 
lective coordinates /3 and 7 pQ. It is also possible to 
parametrize the energy surfaces using the quadrupole 
moments Q 20 and Q 22 related to the /3 and the 7 vari¬ 
ables by /3 = and 7 = tan -1 ^ with Q = 
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Neutron number 


FIG. 9. (Color online) Experimental B(E 2; J —> J — 2) val¬ 
ues in hafnium isotopes (Z=72) (a) and the ratio B 42 = 
B(E2,4i —> 2+)/B(E2, 2^ —> 0^) (b). Values measured in 
this work are marked with a box. Other values are taken from 
the Nuclear Data Sheets I58H62I . 


VQ 20 + Q 22 (see Ref. [33] for more details). In the def¬ 
inition of /3 we use the mean-squared radius (r 2 ) eval¬ 
uated with the corresponding HFB state. Throughout 
this work, the DIM parametrization of the Gogny en¬ 
ergy density functional I .'ill is employed for the effective 
nucleon-nucleon interaction. 

In Fig. [TU] the deformation energy surfaces obtained 
from the Gogny-DIM HFB calculations for the 168 “ 180 Hf 
nuclei, are plotted in terms of the /3 and 7 deforma¬ 
tions. We limit the plot to the range 0.0 < /3 < 0.5 
and 0° < 7 < 60°, because it is the relevant scope for 
our purposes. Energy surfaces for 166 Hf and 180 Hf are 
not shown because they are quite similar in topology to 
the ones of the adjacent nuclei 168 Hf and 178 Hf, respec¬ 
tively. In general, the energy minimum is located around 


TABLE VI. The HFB energy at the minimum (relative to 
the spherical configuration, and denoted by Emin) as well 
as the position of the minimum /32,min are given for both 
parametrizations of the Gogny force and the 170_ 78 Hf nuclei. 



Emin 

D1S 

(MeV) 

DIM 

@2,min 

D1S 

DIM 

TTUjjJ 

-12.214 

-10.944 

0.339 

0.325 

172 Hf 

-13.788 

-12.277 

0.360 

0.332 

174 Hf 

-14.724 

-13.061 

0.353 

0.326 

176 Hf 

-15.068 

-13.385 

0.320 

0.307 

178 Hf 

-14.960 

-13.328 

0.301 

0.288 


/3 = 0.3 on the 7 = 0 axis, being characteristic of an axi¬ 
ally deformed prolate rotor. While any significant change 
in the topology is visible from the microscopic energy sur¬ 
face, the minimum appears to be steeper in both /3 and 
7 directions for heavier Hf isotopes. Although not shown 
here, the results with D1S look quite similar to the DIM 
ones. However, as compared to the DIM results, the en¬ 
ergy minima of the D1S energy surfaces are steeper both 
in (3 and 7 directions than in the DIM case. In Table [Vi] 
we observe that for the nuclei 17 °- 178 Hf the HFB energy 
of the minimum (relative to the energy of the spherical 
configuration (/?, 7)=(0,0) and denoted as E min ) is gen¬ 
erally around 1.2-1.7 MeV smaller in magnitude in the 
DIM case than the D1S one. 

In Fig. [TT] the /3 value at the absolute minimum of the 
microscopic energy surface (denoted hereafter as /32,mm) 
is plotted as a function of neutron number. It exhibits a 
parabolic behavior with its maximum at N = 100 instead 
of the mid-shell value N = 104. The D1S results are 
generally larger than the DIM ones but in both cases 
they show their maximum value at N = 100. However, as 
observed in Table VI the minimum energy E m i n reaches 
its maximum at mid-shell for both the D1S and the DIM 
sets. 

For the boson calculation, we employ the following 
proton-neutron IBM (IBM-2) Hamiltonian: 


Hibm 


e(^d7r d~ rz^) T T zc T • T, (4) 


where the first and the second terms stand for the d- 
boson number operator and the quadrupole-quadrupole 
interactions, respectively. They are given by the expres¬ 
sions n dT = d\-d T and Q T = s\d T + djlsj. +Xt [<4 x d r ]^ 2) 
with r being either 7r (proton) or v (neutron). The third 
term is relevant for rotational bands and it is shown |33| 
to be necessary for the description of deformed nuclei. It 
is given in terms of the total angular momentum operator 
L = Lit + L y with L t = vT0[dt x d T ]^. 

The most general IBM-2 Hamiltonian contains many 
more terms and parameters than the one in Eq. Q. The 
present IBM-2 Hamiltonian is rather simple compared 
with a general Hamiltonian, but contains the minimal 
number of interaction terms relevant for the description 
of low-lying quadrupole states. The Hamiltonian param¬ 
eters are determined by mapping the HFB energy surface 
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/?2 COS7 {3 2 cos 7 P 2 cos 7 


FIG. 10. (Color online) Contour plots of the deformation energy surfaces in terms of the quadrupole deformations /?2 and 7 
for the nuclei 168 - 178 Hf, obtained from the Gogny HFB calculations using the DIM interaction. The color scale varies in steps 
of 0.25 MeV and the contour lines are drawn in steps of 0.5 MeV. The range of the plot is 0.0 ^ /3 ^ 0.5 and 0° 7 ^ 60°. 

The absolute minimum is identified by an open circle. 


onto the IBM one. Some technical details of the proce¬ 
dure are described in Appendix [A} In Table |Vn| we tab¬ 
ulate the IBM-2 parameters determined by the mapping 
procedure for the 166 - 180 nf isotopes. 


TABLE VII. The parameters for the IBM Hamiltonian Hibm 
of Eq. Q, as well as the proportionality coefficient Cp for the 
deformation parameter /3 , obtained from the mapping of the 
HFB to the IBM energy surfaces for 166 - 180 Hf. 


A 

e 

(keV) 

— K 

(keV) 

xlO 3 

xlO 3 

? a 
< ' 

Ce 

eB 

(efm 2 ) 

166 

474 

281 

213 

-813 

-9.95 

3.65 

15.3 

168 

439 

285 

213 

-729 

-11.9 

3.20 

16.3 

170 

451 

292 

241 

-895 

-13.0 

2.96 

15.9 

172 

333 

281 

454 

-878 

-10.5 

2.96 

15.5 

174 

238 

267 

303 

-768 

-6.46 

3.00 

14.3 

176 

166 

242 

237 

-723 

-2.00 

3.30 

12.9 

178 

97.7 

258 

709 

-1086 

-0.952 

3.50 

12.3 

180 

88.0 

273 

359 

-801 

-3.36 

3.60 

13.5 


lated. In particular, the B(E2; J -7 J') value is obtained 
by 

B(E2; J -> J') = ^_\{J'\\f^\\J)\ 2 , (5) 

where J and J' are the total angular momenta of the 
initial and the final states of the transition, respectively. 
The E2^ operator is written as T^ 2 ' 1 = e^Q-* + e^CX,, 
where Q T is the quadrupole operator defined in Eq. (|4j) 
and the same values of X-k,v parameters are used. The 
parameter e T is the boson effective charge for proton and 
neutron. Here we assume that the effective charge is the 
same for protons and neutrons, = e v = es- In most 
of the previous IBM calculations, a fixed value of the 
effective charge es, that is determined by fitting to the 
experimental data for the B{ E2) values, is used for all 
members of the isotopic chain. 

B. Energy levels 


Having all the relevant parameters at hand, the Hamil¬ 
tonian is diagonalized to obtain energies and wave func¬ 
tions of the excited states. For the numerical diagonal- 
ization of the IBM-2 Hamiltonian and the calculations of 
the E2 transition rates, the computer program NPBOS 
f6o| has been used. Using the resultant wave functions, 
the transition probabilities between the states are calcu- 


In order to confirm that the present framework gives a 
reas ona ble description of the energy levels, we display in 
Fig. 12 the experimental level energies for the yrast 2 + , 
4 + , (r, 8 + and 10 + states and the corresponding the¬ 
oretical level energies calculated by the mapped IBM-2. 
Overall, our IBM-2 calculation follows the experimental 
trend satisfactorily, and reproduces the energies for each 


E (MeV) E (MeV) 
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FIG. 11. (Color online) The 02 ,min values, where the HFB 
energy has its absolute minimum are plotted as a function of 
neutron number for the 1G6—ls0 Hf isotopes. Solid and dotted 
curves represent the results from D1S and DIM interactions, 
respectively. 


Experiment IBM2 



Neutron number Neutron number 


FIG. 12. (Color online) Experimental (a) and the theoreti¬ 
cal (DIM parametrization) (b) level energies of the low-lying 
yrast 2 + , 4 + , 6 + , 8 + and 10 + states as functions of neutron 
number. 


nucleus rather well at a quantitative level. However, the 
energy levels have a minimum as a function of neutron 
number at mid-shell N = 104 in the experiment, while 
the minimum is at ./V = 100 in the theory. 

To illustrate this deviation of the 2+ excitation energy, 
we show in Fig. [13] the moments of inertia (denoted as 
MOI) of the rotational band, obtained from the rotor 
formula E(J) oc J(J + 1) (with J = 0 + , 2+, 4 + ,...) and 
the 2 + energies of the IBM result (denoted as “E(2f) 
IBM” in the figure) and the experiment (“E(2+) Exp.”). 
The moment of inertia, which is equal to 3/2£(2+), of the 


FIG. 13. (Color online) The moments of inertia (MOIs) of the 
ground state rotational band of the 166 - 180 jjf nuc lei obtained 
in different approaches are plotted as a function of neutron 
number. The MOIs are computed using the cranking calcu¬ 
lations of the IBM in the coherent state “IBM-cranking” and 
the self-consistent HFB cranking method that leads to the 
Thouless-Valatin MOI both with the DIM “Gogny-HFB TV 
(DIM)” and the D1S “Gogny-HFB TV (D1S)” parametriza- 
tions. Also shown are the MOI extracted from the 2^ excita¬ 
tion energies of the IBM-2 “E(2f) IBM” and the experiment 
“E( 2j") Exp.” using the rigid rotor formula. 


IBM is maximal at N = 100, while the experimental MOI 
at IV = 104 but changes much less than the theoretical 
value with neutron number. 


The reason for the discrepancy in the systematics of 
the MOI between the IBM result and the experiment can 
be due to the inclusion of the L ■ L term in the IBM-2 
Hamiltonian. To shed more light upon this point, we plot 
in the same figure the MOI obtained from the cranking 
calculation in the Gogny HFB model using the Thouless- 
Valatin (TV) formula [M] (denoted as “Gogny-HFB TV” 
in the figure for the two different parametrizations D1S 
and DIM), and the one calculated in the IBM without 
the L ■ L term (denoted as “IBM-cranking”). As de¬ 
scribed in Appendix |aJ the coefficient of the L ■ L term 
is determined so that the “IBM-cranking” MOI becomes 
identical to the “Gogny-HFB TV (DIM)” one. For that 
reason, we do not plot the IBM cranking MOI with the 


L ■ L term. In Fig. 


13 


one should see that the “IBM- 


cranking” MOI is peaked at N = 104 or 106, consistent 
with the experimental MOI while the “Gogny-HFB TV 
(DIM)” MOI is peaked at N = 100 (the same applies 
for the D1S values). For the nuclei with N < 102, the 
difference between the “IBM-cranking” and the “Gogny- 
HFB TV (DIM)” MOIs is rather large as compared to 
the one for N > 104. Thus, the lowering of the energy 
due to the inclusion of the L ■ L is much more significant 
in the N < 102 Hf nuclei than in the N > 104 ones. 
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Consequently, the maximum point of the “E(2f) IBM” 
MOI appears at N = 100 due to the inclusion of the L ■ L 
term. This correlates with the evolution of the derived 
k' parameter, presented in Table VII for example, the 
parameter k' for 172 Hf is much larger in magnitude than 
that for 176-178 Hf. We also note that the same conclusion 
would be extracted from the D1S parametrization as it 
predicts the same systematics as the DIM set. 

It is certainly worth noting that, contrary to the empir¬ 
ical trend, the Thouless-Valatin MOI becomes maximal, 
irrespectively of the choice of the parametrizations, not 
at the mid-shell N = 104 and that this systematics is well 
correlated with the evolution of the /? 2 ,min in Fig. 11 


On the other hand, it has been shown [65] that the 
moment of inertia of the ground-state rotational band 
calculated in the cranking model is rather sensitive to 
the details of the pairing interaction. Therefore, it can 
be of interest that a small modification to the relevant 
channel in the D1S and DIM functionals could lead to a 
substantial improvement of the agreement with the ex¬ 
perimental data. 

Finally, we also note that the calculation generally un¬ 
derestimates the experimental level energies, in particu¬ 
lar, for the nucleus 172 Hf. The reason would be that the 
effect of including the L ■ L term is rather significant, as 
the TV MOI from the Gogny-HFB calculation overesti¬ 
mates the experimental one considerably (see Fig. 131. 


C. E2 transitions 


Let us now turn our attention to the calculation of the 
B( E2) values. In Fig. [l4|a) the theoretical U(E2; J —»• 
J — 2) values (J = 2,4,..., 10) calculated with an ef¬ 
fective charge eg = 0.123 eb are shown. The effective 
charge was fixed for all the considered Hf nuclei as to 
reproduce the experimental B{ E2;2^ —> 0+) value of 
182 W.u in 1,6 Hf. As anticipated from previous IBM 
fitting calculations for W and Os isotopes [5], the the¬ 
oretical -B(E2; J —► J — 2) values with a fixed eg value 
become maximal at mid-shell N = 104, which disagrees 
with the experimental B( E2) systematics showing a peak 
at N = 100 or 102. 

The failure can be partly attributed to the use of a fixed 
effective charge eg. Then, we consider eg to be mass 
dependent and determine its value for each individual 
nucleus. Specifically, we propose to derive eg so that it 
follows the systematics of the Gogny-HFB /? 2 ,min values, 
mainly because this quantity reaches its maximum value 


not at mid-shell but at N = 100, as seen from Fig. 11 


A possible option to take this effect into account can be 
to associate the transition quadrupole moment Q t , corre¬ 
sponding to the 2+ —► 0^ E2 transition matrix element, 
to the Gogny-HFB fa,min value. In general, Q t for the 
transition from the state with spin J to J' is written as 




, 16ttH(E2; J ^ J 1 ) 

~ 5 ~ (J200|J'0) 2 ’ 


( 6 ) 



94 96 98 100 102 104 106 108 

Neutron number 


FIG. 14. (Color online) The B( E2) values for hafnium (Z= 72) 
isotopes calculated with two different choices of effective bo¬ 
son charge eB: (a) eB = 0.123eb fixed for all isotopes; (b) 
es different from nucleus to nucleus, each determined by ad¬ 
justing the transition quadrupole moment Qt(2f — > Oj 7 ) from 
the IBM-2 calculation to fit /?2,mm of the Gogny-HFB DIM 
calculation. The boson effective charges es used for the latter 
calculation are tabulated in Table Em 


with (J200|J'0) being the Clebsch-Gordan coefficient. 
The quantity Qt{J —> J 1 ) is transformed into the in¬ 
trinsic deformation parameter /?t(J —> J') 

A(J ^ J,) = J,) - (7) 

The fit for the 2f —$■ 0] 1 " E2 transition, fit fit y 

is eventually made equal, for each individual nucleus, to 
the Gogny-HFB mean-field fi? min in order to obtain the 
effective charge eg. 

In Fig. [Td])}.) j we plot the resulting B(E2; J — ► J — 2) 
values with the effective charge determined in this way. 
Although the H(E2; J —► J — 2) values in panel (b) are 
generally larger in magnitude than those in panel (a), the 
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systematic trend is more consistent with the experimen¬ 
tal one, as it becomes maximal at N = 100, not at the 
mid-shell nucleus 1,6 Hf. Again, one notices that the re¬ 
sult is well correlated with the Gogny-HFB (3 2 ,min value 
(Fig. Ill, the Thouless-Valatin MOI (Fig. 13), and the 
variation of the extracted ee value as a function of N 


(Tab. VII). 


The observed experimental B{ E2) systematics in the 
neutron-deficient Hf isotopes, indicating the maximum 
collectivity not at the middle of the major shell, implies 
that for the analysis of realistic nuclei a simple model 
may not be sufficient and that certain microscopic effects 
have to be taken into account. In particular, the mass- 
dependence of e-Q implies that the effect of the core polar¬ 
ization may become non-negligible. In the present IBM 
framework, the polarization effect cannot be included ex¬ 
plicitly, but it could be somehow taken into account by 
absorbing it in the variation of the boson charge. 

Alternatively, the mass dependence of the boson ef¬ 
fective charge could be explained by the renormalization 
effect of the g boson [66]. In the sd IBM system, the g- 
boson effect could show up as higher-order terms in the 
quadrupole operator [66] . In that sense, the form of the 
quadrupole operator used in the present work in Eq. Q, 
which is a one-body operator, can be extended to include 
the higher-order terms. 

We also comment on the dependence of the IBM re¬ 
sults on the choice of the EDF parametrization. For the 
energies of the ground-state rotational band, there is a 
certain quantitative difference between the D1S and the 
DIM results, but the overall tendency at the qualitative 
level is expected to be quite similar, because the deforma¬ 
tion energy surface, fi 2 ,min and Thouless-Valatin MOI for 
both parametrizations have been shown to exhibit similar 
features as a function of neutron number. Likewise, for 
the B( E2) systematics, as both the D1S and DIM sets 
provide a similar trend of the f3 2 , m in value (cf. Fig. Ill, 


if the effective charges are determined in the same way, 
one would obtain the results qualitatively similar to the 
ones from the DIM interactions (cf. Figs. 14 and 15). 

For completeness, in Fig. 

£(E2;4+ -> 2+) and B(E2;2+ 


15] we compare the 
0^) experimental val¬ 


ues with the results of the calculations with two different 
choices of the effective charge: (a) e-Q = 0.123 eb fixed 
for all nuclei, and (b) es dependent on mass. In Fig. [l6] 
we compare the experimental and theoretical (with the 
use of mass-dependent effective charge) B( E2;6j l " —► 4+) 
value (a) and the B42 ratio (b). It is rather evident from 
Figs. [l5](a,b) that the calculations with the fixed boson 
effective charge do not reproduce the experimental trend 
while a nice agreement between theory and calculation is 
obtained if one chooses the effective charge for each nu¬ 
cleus as to follow the variation of the /3 2ym in value with N. 
Overall, a reasonable description of the data is obtained 
both qualitatively and quantitatively. 

the calculated B( E2;4/~ 


In Fig. 
B(E2;2+ 


15 


b), 


2J“) and 


O4 ) values overestimate the experimental 


c\j 

LU 

00 
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FIG. 15. (Color online) Comparison of experimental and 
theoretical reduced transition strengths in hafnium isotopes 
(Z=72). (a) Calculation from Fig. [hf](a) (en = 0.123eb). (b) 
(b) (eB mass dependent). 


Calculation from Fig. 


14 


due to the presence of a non-zero hexadecapole deforma¬ 
tion /3 4. A non-zero /I4 has an influence on the transition 
quadrupole moment Qt , and therefore on the extracted 
effective charge. Likely, the inclusion of /I4 will improve 
agreement with experiment for these nuclei. In fact, 
the previous IBM-1 calculation indicated m that the 
hexadecapole degree of freedom may be non-negligible, 
though in a different context of some ground-state prop¬ 
erties. 

For the 64“ —► 4+ E2 transition presented in Fig. 


16 


values for the nuclei around N = 100. This might be 


_ a), 

the theoretical value shows a parabolic behavior. The 
experimental one does not show this trend, but also has a 
large uncertainty. In Fig.[l6|]b), we observe that both the 
experimental and theoretical B42 ratio are close to each 
other, as well as to the SU(3) limit of the IBM (=10/7) 
i! 12] , showing that the considered nuclei are well described 
as being good rotors. 

In comparison to the other EDF-based calculation 
available in the literature, the five-dimensional collective 
Bamiltonian approach based on the Gogny-DlS parame¬ 
ter set has predicted the maximum B(E2;2]/ —► Of) value 
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FIG. 16. (Color online) Comparison of experimental and 
theoretical reduced transition strengths in hafnium isotopes 
(Z=72). (a) Calculation from Fig. [14] (b) (ee mass depen¬ 
dent). (b) B 42 as defined in the text. The IBM SU(3) limit 
of B 42 = 10/7 is indicated as a gray line. 


at 102 f57j. It is one unit different from our result, but 
in that calculation the difference in the _B(E2;2^ —► 0+) 
values between 1 ' 2 Hf and 1,4 Hf is negligible. 


nuclei 172 . 1 74,i76jjj: were remeasured. The results are in 
agreement with those from Coulomb excitation measure¬ 
ments and other methods. The deviation with respect 
to results of 77 fast timing measurements from the late 
1960s can possibly be explained by the improved energy 
resolution of the LaBr detectors and the use of the Or¬ 
ange spectrometer in the current work, both of which 
allow for a better treatment of background and contam¬ 
inations. Furthermore the half lives of the 4+ and 6+ 
state in the nuclei 172 > 174 , 17 6jjf were measured for the 
first time. It was possible to deduce an upper limit for 
the half life of the 8^ state in the same nuclei. 


The evolution of B(E2 ) transition strength in the con¬ 
sidered even-even hafnium isotopes is smooth, as was ex¬ 
pected. The maximum of the B(E 2; 2(j*" —> 0j*~) value 
is found not at mid-shell N = 104, but at lower neu¬ 
tron number, which in this case has turned out to be 
N = 102. A simple model cannot explain this system- 
atics, only predicting the maximum collectivity at the 
mid-shell. One should rather rely on a more realistic or 
microscopic model for complex nuclei. Our spectroscopic 
calculation, performed within the scope of this work, has 
reproduced the experimental trend of the B( E2) very 
nicely (see Fig. [Mjc)), if the boson effective charge is de¬ 
termined as to follow the prediction by the microscopic 
EDF calculation. In the present calculation, any specific 
adjustment to the data has not been invoked, but the re¬ 
sult depends only on the EDF parametrization and on the 
mapping procedure. The results of microscopic EDF cal¬ 
culations on the quadrupole deformation /J 2 min (Fig- H) 
and the cranking moment of inertia (Fig. 13), the corre¬ 
sponding IBM energy levels (Fig. 12) and E2 transition 
rates (Figs. 14 and 15) are consistent and correlated with 
each other very well in systematics. 


On the other hand, the same behavior of the B( E2) 
transition strength has also been observed in neighbor¬ 
ing isotopic chains like erbium, ytterbium and tungsten. 
In that sense, it would be a very interesting subject of a 
future study to analyze the spectroscopy of these neigh¬ 
boring nuclei in a more systematic manner. 


V. CONCLUSION 


The experimental technique of fast timing using a com¬ 
bination of LaBr detectors and a conversion electron 
spectrometer, is very well suited to measure yrast B(E2) 
values in deformed nuclei. It was possible to measure the 
half-lives of ground state band excitations of a nucleus 
from the 2+ to the 6+ state in one single experiment 
with a relatively short amount of measuring time. E.g. 
48 hours of net measuring time in the case of 172 Hf. The 
GCD method for extracting the shorter half-lives gave 
consistent results for different decay-feeder energy com¬ 
binations. The detection limit for picosecond half-lives, 
given a good peak-to-background ratio, was around 5 ps. 

In the present work the half life of the 2+ state in the 
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Appendix A: Procedure to extract parameters for 
the IBM Hamiltonian 


The IBM-2 Hamiltonian of Eq. Q contains 5 free pa¬ 
rameters (e, k, Xv and k') to be determined. The 
bosonic energy surface £tbm(/3,7) is given by an analyt¬ 
ical expression, and is derived by taking the expectation 
value of the IBM-2 Hamiltonian in Eq. Q with the boson 
coherent state |<^(/3, 7 )) [30] : 

|^ l7 )) = H T -^(At)^|0) (Al) 


where N T and |0) represent the number of proton or neu¬ 
tron bosons and inert core, respectively. 


\\ = s* + -^=/3 T sin7 T (d^ +2 + d\_ 2 ) + /3 r cos 7 T d\ 

v’2 


(A2) 


P T is the deformation parameter in the boson system, 
and we assume = p u = ft. In the above equations, 
we distinguish the deformation parameter (3 of the boson 
system from the usual /3 deformation parameter in the 
collective model. The model space of the collective model 
spans the entire nucleus, while only the valence nucleons 
are considered in the IBM. Therefore, the deformation 
parameter for the IBM system ft is always larger than 
the one in the collective model /3, and one can assume, 
to a good approximation, that j3 oc /? [30]. This involves 
an additional parameter Cp, which is the proportionality 
coefficient for the ft deformation, (3 = Cpfi 
The boson energy surface is given as 


Eibm{P, 7) 


<0(/3, 7 )|g IBM |0(/3,7)) 

(<K/3> 7)I<MA 7 )) 

e'(N„ + N v )P 1 2 3 4 N v N v k 
1+P 2 + (l + /3 2 ) 2 


4/3 — 4W -(Xn + Xv)p cos 37 


+ ~X-kXvP 


(A3) 


with e' = e — 6k' and /3 = Cp/3. 

The four parameters e'(= e — 6k'), k, x™ an< 4 Xv plus 
the additional coefficient Cp are determined by adjust¬ 
ing the bosonic energy surface 2?ibm(/3, 7) so that it re¬ 
produces the topology of the microscopic energy surface 
-Ehfb(/3,7) in the neighborhood of the absolute mini¬ 
mum. This procedure reduces to the fitting of the bosonic 
energy surface in Eq. (A3) to the Gogny-DIM energy 


surface. For the fit, we utilize the technique using the 
wavelet transform [31]. 

On the other hand, the L ■ L term does not contribute 


to the energy surface of Eq. (A3) but takes the same ana¬ 
lytical expression as the first term in Eq. ([I]). Therefore, 
once the above five parameters are obtained, the coef¬ 
ficient k' should be determined in a different way from 
the other five parameters. To do this, we take the pro¬ 
cedure of Ref. [33] : the cranking moment of inertia is 
compared between fermion and boson systems. We then 
calculate the MOI for the 2f state by Thouless-Valatin 
(TV) formula [53] : 


Itv = 


J(J+ 1) 
2 £L 


(A4) 


E-y stands for the 2]*~ excitation energy obtained from the 
self-consistent cranking calculation with the constraint 
(J x ) = y/J(J + 1), where J x represents the x component 
of the angular momentum operator. 

The equivalent quantity is derived for the IBM in the 
coherent state |<^>(/3,7)), using the cranking formula of 
Schaaser and Brink [53]: 


7 _ lim 1 <0(1,7) 1^1003,7)) / A c\ 

IBM moow (0(/3,7)|0(/3,7)) 

with u> being the cranking frequency. 

With the parameters e'(= e — 6k'), k, X-k- Xv and Cp, 
fixed from the energy-surface fit, the IBM moment of 
inertia in Eq. (A5) contains only one parameter k'. The 


k! value is determined so that Iibm calculated at the 
equilibrium point, where the energy surface is minimal, 
is equal to the Itv value at the corresponding energy 
minimum. 
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